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ABSTRACT 

We present gas and total mass profiles for 13 low-redshift, relaxed clusters spanning a temperature range 
0.7-9 keV, derived from all available Chandra data of sufficient quality. In all clusters, gas temperature profiles 
are measured to large radii (Vikhlinin et al.) so that direct hydrostatic mass estimates are possible to nearly 
rsoo or beyond. The gas density was accurately traced to larger radii; its profile is not described well by a 
beta-model, showing continuous steepening with radius. The derived p tot profiles and their scaling with mass 
generally follow the Navarro-Frenk- White model with concentration expected for dark matter halos in ACDM 
cosmology. However, in three cool clusters, we detect a central mass component in excess of the NFW profile, 
apparently associated with their cD galaxies. In the inner region (r < 0.1 rsoo), the gas density and temperature 
profiles exhibit significant scatter and trends with mass, but they become nearly self-similar at larger radii. 
Correspondingly, we find that the slope of the mass-temperature relation for these relaxed clusters is in good 
agreement with the simple self-similar behavior, M500 T 01 where a = (1.5 - 1.6) + 0.1, if the gas temperatures 
are measured excluding the central cool cores. The normalization of this M - T relation is significantly, by 
~ 30%, higher than most previous X-ray determinations. We derive accurate gas mass fraction profiles, which 
show increase both with radius and cluster mass. The enclosed / gas profiles within r25oo — 0.4 rsoo have not 
yet reached any asymptotic value and are still far (by a factor of 1 .5 - 2) from the Universal baryon fraction 
according to the CMB observations. The / gas trends become weaker and its values closer to Universal at larger 
radii, in particular, in spherical shells r25oo < r < r^oo. 

Subject headings: cosmology: dark matter — cosmology: observations — X-rays: galaxies: clusters 



1. INTRODUCTION 

Observations of galaxy clusters offer a nu mber of well- 
established cosmological tests (see Voit 2005 for a recent re- 
view). Many of these tests rely on the paradigm in which 
clusters are composed mostly of collisionless cold dark matter 
(CDM), and virialized objects form from scale-free, Gaussian 
initial density perturbations. Numerical simulations of cluster 
formation in CDM cosmology are used to calibrate essential 
theoretical ingredients for cosmological tests, such as the de- 
tailed shape of mas s function models JSheth & Tormer iFl999t 
iJenkins et alJ200lT) or the average baryon bias within clusters. 
The CDM paradigm and numerical simulations make clear 
predictions for the structure of clusters, for example that the 
shoul d have a universal d ensity profile ( Dubinski & Carlber: 
1 199 It iNavarro et ail 1 19961) . that their observable properties 
should exhibit scaling relations, and so on. Confrontation of 
these predictions with the results of high-quality observations 
is a necessary consistency check. Any significant disagree- 
ment, beyond that attributable to variations of an underlying 
cosmology, indicates either that the theoretical models are not 
sufficiently accurate (e.g., they do not include important non- 
gravitational processes) or that there are significant hidden bi- 
ases in the observational studies. Both possibilities are red 
flags for the application of cluster-based cosmological tests in 
the present "era of precision cosmology". 

The above underscores the need for high-quality observa- 
tional studies of representative cluster samples. Of partic- 
ular importance are measurements of the distribution of the 
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dark matter and the dominant baryonic component, the hot 
intracluster medium (ICM). The mass distribution in dynami- 
cally relaxed clusters can be reconstructed using several ap- 
proaches, of which the X-ray method is one of the most 
widely used (a recent review of the mass determination tech- 
niques can be found, e.g., in IVoitl 120051) . X-ray telescopes 
directly map the distribution of the ICM. The ICM in relaxed 
clusters should be close to hydrostatic equilibrium and then 
the spatially-resolved X-ray spectral data ca n be used to de- 
rive the total m ass as a function of radius (e. g jMathewsl 1 9781 
ISarazinll988T) . 

Using X-ray cluster observations for cosmological applica- 
tions has a long history. Early determinations of the amplitude 
of density fluctuations can be found in iFrenk et alJ ( 1990), 
iHenrv & Arnaudl fi99l. iLiiiel i 19921) . and IWhiteetalJ 



( 1993a). Henrv & Arnaucj also used the shape of the cluster 
temperature function to con strain the slope of the p erturbation 
spectrum on cluster scales. Oukbir & Blanchard ( 1992) pro- 
posed that £2,„ can be constrained by evolution of the tempera- 
tu re funct i on, an d thi s test was applied to z ~ 0.3-0.4 clusters 
bv lHenrvl dl997l) andlEke et alJ (119981) . Along a different line 
of argument. IWhite etaiT ( [T993bh obtained the first determi- 
nation of Q,„, assuming that the baryon fraction in clusters ap- 
proximates the cosmic mean. Attempts to use the gas fraction 
as a d ista nce indicator were mad e in papers by Rine s et all 
( 1999) and Ettori & Fabian ( 1999). The common limitation of 
these early studies is the rather uncertain estimates of the clus- 
ter total mass. Accurate cluster mass measurements at large 
radii are challenging with any technique. The X-ray method, 
for example, requires that the ICM temperature be measured 
locally in the outer regions; the hydrostatic mass estimate is 
only as accurate as T and dT/dr at radius r. 

Spatially resolved cluster temperature measurements first 
became possible with the launch of ASCA and Beppo-SAX 



2 



VIKHLININ ET AL. 



TABLE 1 
Cluster Sample 



Cluster Z r mm 

" r dct b ROSAT* 



A133 0.0569 40 1100 + 

A262 0.0162 10 450 + 

A383 0.1883 25 800 

A478 0.0881 30 2000 + 

A907 0.1603 40 1300 

A1413 0.1429 20 1800 

A1795 0.0622 40 1500 + 

A1991 0.0592 10 1000 + 

A2029 0.0779 20 2250 + 

A2390 0.2302 80 2500 

RXJ1159+5531 0.0810 10 600 

MKW4 0.0199 5 550 + 

USGCS152 0.0153 20 300 



11 — Inner boundary (kpc) of the radial range used for the tem- 
perature profile fit flBl. 

b — The radius (kpc) where X-ray brightness is detected at 
> 3<x, or the outer boundary of the Chandra field of view. 

c — Those clusters for which we use also ROSAT PSPC surface 
brightness measurements. 

(e.g.. iMarkevitch et alJ fl998t iDe Grandi & Molendil 120021 
Truly accurate temperature profiles are now provided by 
Chandra and XMM-Newton. Chandra results on the mass 
dis tribution in t he inn e rmost cluste r regio n s were re ported 
bv iDavid et al] J2001I). iLewis et alJ J2003I) . iBuote & Lewisl 
(2004), and Arabadjis et al. (2004), and within ~ 1/4 of r2oo 
in a series of papers by lAllen et all (|2004, and references 
therein). The XMM-Newton mass measurements at large 
radii (~ 0.5 r ono) in 10 low-redshift cl usters we re recently 
publis hed by IPointecout eau et all J2005I) and lArnaud et alJ 
<200i>. 

In this paper, we present the mass measurements in a sam- 
ple of 13 low-redshift clusters whos e temperature profiles 
were derived in JVikhlinin et alJl200 5. Paper I hereafter; Ta- 
ble 0. These clusters were observed with sufficiently long 
Chandra exposures that the temperature profiles can be mea- 
sured to 0.75 rsoo (~ 0.5 r2oo) in all objects and in 5 cases, 
can be extended outside r^oo- All these objects have a very 
regular X-ray morphology and show only weak signs of dy- 
namical activity, if any. Even though the present sample is not 
a statistically complete snapshot of the cluster population, it 
represents an essential step towards reliable measurements of 
the cluster properties to a large fraction of the virial radius. 

The Paper is organized as follows. Our approach to 3- 
dimensional modeling of the observed projected quantities is 
described in §[3j Results for individual clusters are presented 
in §|4]and self-similarity of their temperature and density pro- 
files is discussed in §|5] and |6] Our data lead to accurate de- 
termination of the M - T relation for relaxed clusters (§0. 
In §|8j we discuss the observed systematic variations of ICM 
mass fraction with radius and cluster mass. 

To compute all distance-dependent quantities, we assume 
Qm = 0.3, A = 0.7, h = 0.71. Uncertainties are quoted 
at 68% CL. Cluster masses are determined at radii rsoo and 
r25oo, corresponding to overdensities 500 and 2500 relative to 
the critical density at the cluster redshift. 

2. X-RAY DATA ANALYSIS 

The main observational ingredients for the present analy- 
sis are radial profiles of the projected temperature and X-ray 
surface brightness. We refer the reader to Paper I for an exten- 
sive description of all technical aspects of the Chandra data 




Fig. 1. — Detected sources in the Chandra ACIS-I field of A1413. Yellow 
circles mark point sources. The only detectable extended X-ray source (other 
than A1413) is marked by the red circle. ACIS-I field of view is 16' X 16'. 



reduction and spectral analysis and discuss here only the X- 
ray surface brightness profile measurements. 

First, we detected and masked out small-scale X-ray 
sources detectable in either soft (0.7-2 keV) or hard (2-7 ke V) 
energy bands. Detection was perfor med using the wavele t 
decomposition algorithm described in Vikhlinin et al. ( 1998). 
Detection thresholds were chosen to allow 1-2 false sources 
per field of view. The results were hand-checked for each 
cluster and exclusion radii adjusted if needed. Source de- 
tection for one of our clusters (Fig.^ is illustrated in Fig.^ 
Note that in addition to point sources, we detected and ex- 
cluded small-scale extended X-ray sources (a typical exam- 
ple of a source of this type is shown by the red circle in 
Fig-E- This should remove at least the brightest of the cold 
gas clumps associated with groups and individual galaxies 
which are p resent in the cluster volume jMotl et all l2004t 
Nag ai et alJl2003l iDolag et alJ l2004b) and could bias mea- 
surements of the global cluster parameters jRasia et alJ2005l) . 
The typical limiting flux for detection of compact extended 
sources is ~ 3 x 10" 15 erg s" 1 crrT 2 in the 0.5-2 keV band; 
this corresponds to a luminosity of ~ 1 .5 x 10 42 erg s for the 
median redshift of our sample, z = 0.06. 

The surface brightness profiles are measured in the 0.7- 
2 keV energy band, which provides an optimal ratio of the 
cluster and background flux in Chandra data. The blank-field 
background is subtracted from the cluster images, and the re- 
sult is flat-fielded using exposure maps that include correc- 
tions for CCD gaps and bad pixel, but do not include any spa- 
tial variations of the effective area. We then subtract any small 
uniform component corresponding to adjustments to the soft 
X-ray foreground that may be required (see Paper I for details 
of background modeling). This correction is done separately 
for the back-illuminated (BI) and front-illuminated (FI) CCDs 
because they have very different low-energy effective area. 

We then extracted the surface brightness profiles in narrow 
concentric annuli (r out /r; n = 1.05) centered on the cluster 
X-ray peak and computed the area-averaged Chandra effec- 
tive area for each annulus (see Paper I for details on calcu- 
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lating the effective area). Using the effective area and ob- 
served projected temperature and metallicity as a function of 
radius, we converted Chandra count rate in the 0.7-2 keV 
band into integrated emission measure, EM - Jn e n p dV, 
within the cylindrical shell. We have verified that these "phys- 
ical" cluster brightness profiles derived from BI and FI CCDs 
are always in excellent agreement in the overlapping radial 
range. The profiles from both CCD sets and different point- 
ings were combined using the statistically optimal weighting. 
Poisson uncertainties of the original X-ray data were propa- 
gated throughout this procedure. 

For lower-redshift clusters in our sample, the statistical 
accuracy of the surface brightness at large radii is limited 
mostly by the Chandra field of view. The analysis in such 
cases can benefit from also using data from the ROSAT 
PSPC pointed observations. We used flat-fielded ROSAT im- 
ages in t he 0.7-2 keV band p roduced using S. Snowden's 
software JSnowdenetalJll994l) and reduced as described in 
IVikhhnin et all J1999I) . ROSAT surface brightness profiles 
were converted to projected emission measure and were used 
as a second independent dataset in modeling the gas density 
distribution. We did not use the ROSAT data in the central 2' 
from the cluster center, because this region can be affected by 
the PSPC angular resolution (~ 25" FWHM). At larger radii, 
we always find excellent agreement between Chandra and 
ROSAT PSPC surface brightness data (an example is shown 
in Fig.|2}. Below, we used the combined Chandra and ROSAT 
analysis for those clusters with sufficiently deep PSPC ex- 
posures — A133, A262, A478, A1795, A1991, A2029, and 
MKW4. 

3. MODELING OF TEMPERATURE AND SURFACE 
BRIGHTNESS PROFILES 

We model the observed X-ray surface brightness and pro- 
jected temperature profiles using the following general ap- 
proach. The 3D profiles of gas density and temperature are 
represented with analytic functions which are smooth, but 
have freedom to describe a wide range of the possible profiles. 
The models are projected along the line of sight and fit to the 
data. The best fit 3D model is used to derive all interesting 
cluster parameters, such as the total gravitating mass. Mea- 
surement uncertainties for all quantities are estimated using 
Monte-Carlo simulations in which random statistical errors 
are added to the data and the full analysis is repeated using 
these simulated data as an input. 

Reliability of this modeling approach was tested by apply- 
ing our Chandra data analysis procedures to "observations" of 
clusters from high-resolution numerical simulations (Nagai et 
al., in preparation). This analysis demonstrated that the three- 
dimensional gas density and temperature profiles of relaxed 
clusters are reconstructed within a few percent. 

3.1. Gas Density Model 

The analytic expression we use for the 3D gas density dis- 
tributio n is obtained by m odifying the traditionally used f3- 
model (Cavaliere & Fusco-Femiano 1978). Modifications are 
designed to represent essential features of the observed X-ray 
surface brightness profiles. Gas density in the centers of re- 
laxed clusters, such as ours, usually has a power law-type cusp 
instead of a flat core and so the first modification is 

< _^ < W^)'* m 

(1 + r 2 lr 2 fP (l+r 2 /r 2 c )W-*/ 2 ' 1 ' 

This model was also used by, e.g.. lPointecouteau et al.N2004l) . 
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Fig. 2. — Observed projected emissivity profile for A133. Chandra and 
ROSAT PSPC data are shown in red and green, respectively. Solid line shows 
the best fit to the 3D gas density model (eq. |3)). The slope of the emissivity 
profile steepens by 1 at radius rt, r = r,(e - l) _1 ' r . 

At large radii, the observed X-ray brightness profiles of- 
ten steepen at r > 0.3 r? nn relative to the powe r law extrapo- 
lated from smaller radii ( Vik hlinin et aljfl999h . This change 
of slope can be modeled as 



1 



(1 +r 2 /r 2 F 



(i + r 2 /r 2 ) 3 P (i + fjr s yyiy ' 



(2) 



where the additional term describes a change of slope by s 
near the radius r s and the parameter y controls the width 
of the transition region. Finally, we add a second /3-model 
component with small core-radius to increase modeling free- 
dom near the cluster centers. The complete expression for the 
emission measure profile is 



(r/r c y 



1 



(1 + r 2 /r 2 ) 3 /?- a / 2 (1 + rvlr s vyly 
n 2 

"02 



(3) 



(1 + r 2 /r 2 2 ) 3 A ' 



All our clusters can be fit adequately by this model with a 
fixed y — 3. All other parameters were free. The only con- 
straint we used to exclude unphysically sharp density breaks 
was s < 5. The analytic model has great freedom and 
can fit independently the inner and outer cluster regions. This 
is important for avoiding biases in the mass measurements 
at large radii and also for more realistic uncertainty estimates. 
The best-fit model for the observed surface brightness in Al 33 
is shown by a solid line in Fig. [2] 

Parameters in eq. are strongly correlated and therefore 
their individual values are degenerate. This is not a problem 
because our goal is only to find a smooth analytic expression 
for the gas density which is consistent with the observed X- 
ray surface brightness throughout the radial range of interest. 

3.2. Temperature Profile Model 

Several previous studies used a polytropic law, T(r) cc 
[Pgas ( r) ] y -1 , to model n o n-con st ant cluster temperature 
pro files JMarkevifch et alJ 119991 iFinoguenov et al] 2001; 
iPratt & Am aud 2002]). We take a different approach because 
the polytropic model is in fact a poor approximation of the 
temperature profiles at large radii jMarkevitch et alJ Il998t 
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TABLE 2 

Best Fit Parameters for Gas Density Profiles 



Cluster 


relet " 


«o 


r c 


r, 


(X 


P 


g 


"02 




Pi 


ROi>AT° 




(kpc) 


10 J cm 3 


(kpc) 


(kpc) 








10 1 cm 1 








A133 


1 100 


2 968 


142 7 


1423 3 


996 


575 


5 000 


276 


33 44 


980 




A262 


450 


3.434 


45.2 


350.8 


1.674 


0.333 


1.806 








+ 


A383 


800 


7.000 


115.2 


422.3 


2.018 


0.583 


0.740 


1.014 


0.08 


1.000 




A478 


2000 


8.169 


177.1 


3148.2 


1.493 


0.715 


5.000 


0.584 


24.00 


1.000 


+ 


A907 


1300 


6.257 


136.9 


1885.1 


1.554 


0.594 


4.986 










A1413 


1800 


5.526 


186.3 


2077.1 


1.217 


0.651 


4.991 










A1795 


1500 


14.993 


72.8 


1030.8 


1.060 


0.545 


3.474 








+ 


A1991 


1000 


9.373 


44.4 


998.2 


1.516 


0.501 


5.000 


0.999 


5.00 


1.165 


+ 


A2029 


2250 


16.469 


80.7 


870.6 


1.131 


0.539 


1.650 


3.741 


5.00 


1.000 


+ 


A2390 


2500 


3.069 


353.6 


1200.0 


1.917 


0.696 


0.240 










RXJ 1159+5531 


600 


0.198 


613.8 


961.5 


1.762 


1.215 


4.939 


0.416 


12.66 


1.000 




MKW4 


550 


0.280 


488.6 


1081.6 


1.628 


1.224 


0.000 


0.189 


11.08 


0.661 


+ 


USGCS152 


300 


17.450 


8.1 


467.5 


2.644 


0.453 


3.280 











Note. — Derived densities and radii scale with the Hubble constant as h 1 ' and h , respectively. 

a — The radius (kpc) where X-ray brightness is detected at > 3o\ or the outer boundary of the Chandra field of view, whichever is smaller. 
b — Those clusters for which we use also ROSAT PSPC surface brightness measurements. 



iDe Grandi & Molendil l2002). which is apparent in our more 
accurate Chandra measurements. 

The projected temperature profiles for all our clusters show 
similar behavior (see Fig. 15 in Paper I and Fig.l3Ul5lbelow). 
The temperature has a broad peak near 0. 1-0.2 of r2oo and de- 
creases at larger radii in a manner consistent with t he results of 
earlie r ASCA and Beppo-SAX obser vations JMarkevitch et al.l 
119981 IDe Grandi & Molendill2002T) . reaching approximately 
50% of the peak value near 0.5 r2oo- There is also a tempera- 
ture decline towards the cluster center probably because of the 
presence of radiative cooling. We construct an analytic model 
for the temperature profile in 3D, in such a way that it can 
describe these general features. Outside the central cooling 
region, the temperature profile can be adequately represented 
as a broken power law with a transition region: 

(1 +(r/r t ) b y/ b ' 

The temperature decline in the central region in most clusters 
can be described as 

fcooiM = (x + aWJ-oVOc +1), x = (r/r cool ) fl ™' , (5) 

JAllen. Schmidt & Fabian 1200 lbT) . Our final model for the 3D 
temperature profile is the product of eq. |3] and (5), 

TiD(r) = To x f coo i(r) x t(r). (6) 

Our model has great functional freedom (9 free parameters) 
and can adequately describe almost any type of smooth tem- 
perature distribution in the radial range of interest. 

This model is projected along the line of sight to fit the 
observed projected temperature profile. This projection re- 
quires proper weighting of multiple tem perature components. 
We use the algorithm described in Vikhlinin (2005|) that 
very accurately predicts the single-temperature fit to multi- 
component spectra over a wide range of temperatures. The 
inputs for this algorithm are the 3D profiles of the gas temper- 
ature, density, and metallicity. The only missing ingredient 
in our case is the 3D metallicity profile. We use instead the 
projected metallicity distributions presented in Paper I. The 
difference between the projected and 3D abundance profiles 
leads to very small corrections in the calculation of the project 
temperature that are negligible for our purposes. We also note 
that the 3D temperature fit at large radii is rather insensitive 



to the choice of weighting algorithm. For example, we tested 
the commonly used method of weighting T with the square 
of the gas density and obtained very similar results. The pri- 
mary reason for stability of the r pl0 j calculations is that the 
ICM emissivity is a strongly decreasing function of radius and 
most of the emission observed at a projected distance b comes 
from a narrow radial range near r — b. 

In several cases, the cluster X-ray brightness is detected 
beyond the outer boundary of the Chandra temperature pro- 
file. Typical examples are A262 where the ROSAT PSPC pro- 
file extends to 50', well beyond the Chandra FOV, and A383 
where the X-ray brightness is detectable in the Chandra im- 
age to 1300 kpc, while the temperature profile is sufficiently 
accurate only within the central 750 kpc. Detection of ICM 
emission at large radii sets a lower limit to the temperature 
and thus provides additional information for the T(r) mod- 
eling. We required that the 3D temperature model exceeds 
0.5 keV at r — r^ el , the radius where the X-ray brightness is at 
least 3cr significant (rd et are reported in Table^. 

We excluded from the fit the data within the inner cutoff ra- 
dius, TYnjn (listed in Tabled, which was chosen to exclude the 
central temperature bin (10-20 kpc) because the ICM is likely 
to be multi-phase at these small radii. In A133 and A478, r m j n 
was increased to exclude substructures associated with activ- 
ity of the central AGN. The cutoff radius was increased also in 
A1795, A2390, and USGC S152, because our analytic model 
is a poor fit to the inner temperature profile in these clusters. 
The choice of r^n is unimportant because we are primarily 
interested in the cluster properties at large radii. 

3.3. Uncertainties 

Our analytic models for the gas density and temperature 
profiles have many free parameters and strong intrinsic degen- 
eracies between parameters. Therefore, uncertainty intervals 
for all quantities of interest were obtained from Monte-Carlo 
simulations. The simulated data were realized by scattering 
the observed brightness and temperature profiles according to 
a Gaussian distribution with a dispersion equal to measure- 
ment uncertainties. The surface brightness and temperature 
models were fit to the simulated data and the full analysis re- 
peated. The uncertainty on all quantities of interest were ob- 
tained by analyzing the distribution derived from 1000-4000 
simulated profiles we generated for each cluster. 
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Fig. 3. — Results for A133. (a) — Temperature profile. Observed projected temperatures are shown by crosses. Solid red and blue lines show the best fit 3D 
model and the corresponding projected profile, respectively. Dotted lines indicate the 68% CL uncertainty interval obtained from Monte-Carlo simulations (see 
text). Models are shown in the radial range r m j n -r,j ct (see text), (b) — Mass and density profiles. M(r) increases with radius and p(r) decreases. Red and blue lines 
show results for the total mass and gas mass, respectively, (c) — Gas mass fraction as a function of radius. Lines show the enclosed / eas = M gas (< r)/M tot (< r). 
The local gas fraction (p ga s/Ptot) in the radial range directly covered by the Chandra temperature data is shown by crosses. The vertical line shows the radius rsoo 
derived from the best fit mass model. 
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Fig. 4. — Results for A262. See caption to Fig.l3l 
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3.4. Mass Derivation 

Given 3D models for the gas density and temperature pro- 
files, the total mass within the radius r ca n be esti mated from 
the hydrostatic equilibrium equation (e.g., Sarazin 1988): 

n (dlogp„ d\oeT\ 

M(r) = -3.7 lxlO 13 M r(r) r — ^ + — 2- , (7) 

\ a log r a log r ) 

where T is in units of keV and r is in units of Mpc. Given 
M(r), we can calculate the total matter density profile, p(r) = 
{Anr 2 y l dM/dr. We also compute the total mass at several 
overdensity levels, A, relative to the critical density at the clus- 
ter redshift, by solving equation 

M A (r A ) = A 4/3 nrlpAz), (8) 

(see, e.g.. lWhite!200!l for discussion of different cluster mass 
definitions). 

The ICM particle number density profile is given directly 
b y the analytic fit to the projected emission measure profile 
(S I3.lt and it is easily converted to the gas density. For the cos- 
mic pl asma with the He abundance from Anders & Grevesse 
(1989),p g = 1. 252 m p (n p n e ) 1/2 . 

The total mass derived from eq. Q is a complex combina- 
tion of parameters which define our T(r) and p g (r) models. Its 
uncertainties are b est derived via Monte-Carlo simulations as 
described in § 13.31 We used the best fit T(r) model obtained 
for each realization of the temperature profile to compute all 



quantities of interest (Ma, M g , gas mass fraction, etc.). The 
peak in the obtained distribution corresponds the most proba- 
ble value ("best fit") and the region around the peak contain- 
ing 68% of all realizations is the 68% CL uncertainty interval. 

Our analytic model for T(r) allows very steep gradients. In 
some cases, such profiles are formally consistent with the ob- 
served projected temperatures because projection washes out 
steep gradients. However, large values of dT/dr often lead 
to unphysical mass estimates, for example the profiles with 
p < at some radius. We eliminated this problem in the 
Monte-Carlo simulations by accepting only those realizations 
in which the best-fit T(r) leads to p tot > p gas in the radial range 
covered by the data, r m \ n < r < r^. Finally, we checked 
that the temperature profiles corresponding to the mass uncer- 
tainty interval are all convectively stable, din T/dlnp g < 2/3. 

3.5. Average Temperatures 

We also computed average temperatures for each cluster us- 
ing different weightings of the 3D temperature models. The 
temperatures were averaged in the radial range 70 kpc < r < 
r5oo- The central 70 kpc were excluded because temperatures 
at these radii can be strongly affected by radiative cooling and 
thus not directly related to the depth of the cluster potential 
well. The averages we compute are: 

T m g — weighted with p gas (r). T m „ is needed, e.g., to com- 
pute the integrated SZ (Sunvaev & Zeldovichl ll972l) signal. 
T mg should be more directly related to the cluster mass than 
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Fig. 7. — Results for A907. See caption for Fig.l3l 



X-ray emission-weighted temperatures. 

T^spec — a value that would be derived from the single- 
temperature fit to the total cluster spectrum, excluding the 
central region. r spec is obtained by integrating a combination 
of T(r) and p\. is {r) as described in IVikhlininl d2005l) . r spec is 
one of the primary X-ray observables. 

Measurement uncertainties for the av erage te mper atures are 
derived via Monte-Carlo simulations (§ !3.3l and l3.4> . 

4. RESULTS FOR INDIVIDUAL CLUSTERS 

Best fit parameters for the gas density distribution are listed 
in Table|3] We do not provide measurement uncertainties for 
the individual parameters because they are strongly degener- 



ate. Monte-Carlo simulations show that for each cluster, the 
statistical uncertainties of the gas density are below 9% every- 
where in the radial range 5 kpc - r^. However, extrapolations 
of the gas density far beyond r^ et are unreliable. The main de- 
rived cluster parameters are reported in Table[3] and tempera- 
ture, mass, and gas fraction profiles for individual clusters are 
presented in Fig.l3TH31 

The observed temperature profile and the best-fit models 
are shown in panels (a). The models are shown in the radial 
range r m ; n - r^. Note that our analytic function successfully 
describes a wide range of shapes for the observed temperature 
profiles and the corresponding projected profiles (shown by 
blue lines) are always excellent fits to the data. The 3D models 





(red lines) go above the data in the outer region, and below the 
data in the center — just as expected from projection of the 
temperature gradients along the line of sight. 

Dotted lines show 68% CL uncertainties for the model at 
each radius derived from Monte-Carlo simulations. The un- 
certainty intervals are smaller than the errorbars of the raw 
measurements because the model effectively smoothes the 
data over 3-4- adjacent bins (and therefore uncertainties in the 
neighboring bins are correlated). However, the difference is 
not very large and the derived uncertainties are clearly real- 
istic in the sense that they include a typical range of smooth 
models that could be drawn through the data. 

Panels (b) present the derived density and enclosed mass 



profiles. Results for the total and gas masses are shown by 
red and blue lines, respectively. M(r) increases with radius 
and p(r) decreases. The derived total mass profile is used to 
estimate the radius, rsoo (vertical dotted lines). If rsoo is out- 
side the radial range covered by the X-ray data, we find it 
from extrapolation of our best fit gas density and/or tempera- 
ture models. Large extrapolations are required only for A262, 
RXJ 1159+5531, and USGC S152. In other cases, Chandra 
temperature measurements extend almost to rsoo or beyond. 

Dotted lines show 68% CL uncertainties for the mass and 
density derived from Monte-Carlo simula tion s for a subset of 
physically meaningful realizations (see § !3.4t . Uncertainties 
are shown also for the gas mass but they are very small and 




almost indistinguishable from the best fit model in these plots. 

Gas fraction profiles are presented in panels (c). The en- 
closed gas fraction is shown by red lines. The local gas 
fraction, p gas /ptot, in the radial range directly covered by the 
Chandra temperature data, is shown by crosses. 

Abell 2390 deserves a special note. Its deep Chandra im- 
age reveals large-scale cavities in the X-ray surface brightness 
extending ~ 400 kpc from the center where a shar p break 
in the s urface brightness profile was reported by Allen et al. 
(2001a). The cavities are likely produced by bubbles of radio 
plasma emitted by the central AGN, as observed in several 
other clusters (McNamar a et al.ll2005l) . The ICM in the cen- 
tral region is not spherically symmetric nor expected to be in 



hydrostatic equilibrium. This should result in underestimation 
of the total mass and overestimation of the gas mass. The re- 
sults at small radii for A2390 should be treated with caution. 
There are no detectable structures outside 500 kpc and so the 
results at large radii (e.q., at rsoo) should be reliable. 

5. AVERAGE TEMPERATURE PROFILE 

We noted in Paper I that the temperature profiles for our 
clusters are self-similar when scaled to the same overden- 
sity radius, in good ag reem ent with the earlier studies by 
Mark evitch et alJ lQj98) and [De Grandi & Molendil <2002t . 
We return to this subject here because results for individual 
clusters can be compared more accurately using the recon- 
structed 3 -dimensional temperature profiles and the overden- 




sity radii determined from the mass model rather than esti- 
mated from the average temperatures as in Paper I. 

Figure ^] shows the reconstructed 3-dimensional temper- 
ature profiles, normalized to the gas mass-weighted average 
temperature, T mg , and plotted as a function of r/rsoo- The 
model for each cluster is plotted in the radial range directly 




T < 2.5 keV 
A2390 



Fig. 16. — Scaled 3-dimensional temperature profiles. The mean temper- 
ature is the gas mass-weighted average. The thick yellow line shows an ap- 
proximation to the average profile for T > 2.5 keV clusters (eq. 1^1). 



covered by the Chandra spectral measurements and our 3D 
modeling, from r m m to the outer bin in Fig.l3*h-tl5k. We nor- 
malize the profiles to the gas mass-weighted temperature be- 
cause it is less sensitive than, e.g., T SJ)ec to the properties of 
the central regions where the ICM temperature can be signifi- 
cantly affected by non-gravitational processes. 

At r > 0.05 rgoo, the scaled 3-dimensional temperature for 
all but one of our T > 2.5 keV clusters are within +15% from 
the average profile which can be approximated with eq. Q, 

7\r) = ^ (x/0.045)'- 9 + 0.45 1 

T mg ' (jt/0.045) 1 - 9 + 1 (1 + (x/0.6) 2 ) 45 ' 

where x = r/rsoo- This average model is shown in Figll6lbv 
thick yellow line. The only deviation is A2390, in which the 
inner temperature profile and overall normalization appear to 
be distorted by activity of the central AGN (see above). The 
profiles of the low-temperature clusters, T < 2.5 keV (green 
lines in Fig.ll6t. follow the average profile at r > 0.3 rsoo but 
are significantly different and show large scatter at small radii. 
They show a stronger temperature increase to the center, and 
the central cooler regions seem to be confined to a smaller 
fraction of the virial radius than in the more massive clusters. 

The observed radial temperature variations imply that the 
average cluster temperature cannot be defined uniquely. A 
possible difference between different definitions (spectro- 
scopic, emission-weighted, gas mass-weighted etc.) should 
be kept in mind. Also, the aperture size used for integration of 
the X-ray spectrum is important. Our average 3-dimensional 
profile (5) implies the following approximate relation between 
the peak, spectroscopic average, and gas mass-weighted tem- 
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TABLE 3 

Masses, Average Temperatures, Concentration parameter, Gas Fractions 



Cluster rsoo ^spec T mg C500 M15Q0 M500 /g, 2500 fg, 500 /g, 2500-500 

(kpc) (keV) (keV) (10 14 M o ) (10 14 M o ) 

A133 998±37 4.15 ±0.07 3.68 ± 0.1 1 3.31±0.27 1.24±0.08 3.14±0.36 0.065 ± 0.002 0.084±0.006 0.095^g 

A262 663 ± 19 2.08 ± 0.06 1.92 ±0.09 3.58 ± 0.27 0.34 ± 0.05 ■■■ 0.067 ± 0.003 ■■■ ■■■ 

A383 956 ±33 4.80 ±0.12 4.36 ±0.18 4.36 ± 0.35 1.68 ±0.15 3.10 ±0.32 0.090 ± 0.005 0.122 ±0.007 0.149+^ 

A478 1359±60 7.95 ±0.14 7.34 ± 0.32 3.75 ± 0.28 4.23 ± 0.26 7.83 ± 1.04 0.096 ± 0.004 0.118 ± 0.011 0. 1 30+°,°,^ 

A907 1117 ± 31 5.96±0.08 5.44±0.13 3.56±0.31 2.30±0.16 4.71±0.39 0.090 ± 0.003 0.121 ±0.005 0.145 + ° : ™ 

— u.Olb 

A1413 1339 ±48 7.38 ±0.12 6.76 ± 0.20 3.06 ±0.16 3.08 ±0.19 7.78 ± 0.83 0.092 ± 0.003 0.104 ± 0.007 0.1 12+ggj^ 

A1795 1283 ± 46 6.10 ±0.06 5.52 ±0.11 3.73 ±0.17 2.66 ±0.15 6.57 ± 0.69 0.088 ± 0.003 0.098 ± 0.007 0.099^ggJ^ 

A1991 753 ±39 2.59 ± 0.06 2.23 ± 0.09 4.66 ± 0.32 0.60 ± 0.07 1.28 ±0.20 0.068 ± 0.004 0.099 ± 0.009 0.103+^ 

A2029 1380 ±31 8.46 ± 0.09 7.59 ±0.19 4.13 ±0.24 4.41 ±0.23 8.29 ± 0.79 0.090 ± 0.003 0.121 ±0.007 0.144+^3 

A2390 1448 ±49 8.90 ±0.17 9.35 ± 0.39 1.65 ±0.18 3.50 ± 0.28 10.88 ± 1.05 0.127 ±0.005 0.138 ±0.008 0.141+^];} 

RXJ 1 159+5531 ... 695 ±60 1.80±0.10 1.58±0.09 1.89±0.35 0.31 ±0.04 ■■■ 0.042 ± 0.002 ■■■ ■■■ 

MKW4 635 ±27 1.65 ±0.05 1.40 ±0.05 2.69 ±0.17 0.30 ± 0.03 0.74 ± 0.09 0.044 ± 0.002 0.063 ± 0.005 0.073+° [Jj 4 

USGCS152 ■■■ 0.69 ±0.02 0.59 ± 0.03 ••• 0.07 ± 0.00 ■■■ 0.043 ± 0.001 ■■■ ■■■ 



Note. — The derived quantities scale with the Hubble constant as rjoo K h , M oc h , /g as oc K 



peratures (all measured in the radial range 70 kpc - rsoo): 

r peak : r spec : T mg = 1.21 : 1.11 : 1. (10) 

Comparing our temperature profi les with the c ompilation 
of XMM-Newton results presented in Arnaud et al. ( 2005), we 
note a general agreement at small radii. However, the results 
at large radii seem to be different — the temperature decline 
is generally not present in these XMM-Newton clusters. Dis- 
cussion of this discrepancy is beyond the scope of this paper. 
The arguments for validit y of our m easurements and detailed 
comparison with the Arnaud et al. results for several clusters 
in common can be found in P aper 1 

6. TOTAL AND GAS DENSITY PROFILES 

One of the key theoretical predictions of the hierarchi- 
cal CDM models is the universal density di stribution within 
dark matter halos (iNavarro et alJl996lH997l hereafter NFW). 
Specifically, the shape of the radial density profiles of Cold 
Dark Matter halos is characterized by a gradually changing 
slope a = <ilogp/<ilog r from a « -1 in the inner regions 
to a ~ -3 at large radii (Dubinski & Carlber jlliWll NFW). 
The profiles are characterized by concentration, ca, defined 
as the ratio of the halo virial radius and the scale radius, r s : 
ca = r&/r s . The scale radius is defined as the radius where the 
logarithmic slope of the density profile is a = -2. Concen- 
trations of CDM halos are tightly correl ated with the c harac- 
teristic epoch of object formation (Wechsler et al. 2002J). The 
mean concentration is only a wea kly decreasing function of 
the virial mass, c\ cc M " 1 , JNavarro et al.ll997tlBullock et alJ 
l27MlEkeetani200l . Therefore, within a limited range of 
masses, the mean concentration is approximately constant and 
the density profiles are approximately self-similar. 

Such self-similarity is indeed observed in our data. In 
FigEl we pl°t the derived total density profiles, scaled to 
p c (z), as a function of radius in units of rsoo- The scatter of 
individual clusters around a mean profile is small and consis- 
tent with that found in numerical simulations (see § 16. ll be- 
low). The average total density profile in our clusters agrees 
well with the NFW model with the concentration expected for 
objects of this mass (thick yellow line in Figll7>. 



The p tot profiles for 3 clusters with T < 2.5 keV have cen- 
tral steepening at r < 0.05 rsoo ~ 30 kpc, which is statistically 
significant. Such sharp steepening is not expected in purely 
CDM halos. We associate these components with the stel- 
lar material of the central cD galaxies, which start to domi- 
nate the total mass at small radii. Indeed, our derived masses 
within 30 kpc for these objects, ~ 2 x 1O 12 M (see Fig0 
are similar to the stel lar mass es timates in large central cluster 
galaxies (e.g.. lLin & Mohrl2004l) . 

Self-similarity is also observed in the ICM distribution at 
large radii. The lower set of profiles in Fig.[n]shows the gas 
densities, also scaled to p c {z). The scaled gas density near rsoo 
is within +15% of the mean for most clusters. Some of this 
scatter is caused by uncertainties in the total mass estimates 
— a typical uncertainty of 3-5% in rsoo (Table |3j translates 
into 7-12% scatter in the scaled densities for a typical slope 
of PgasM near rsoo- The scatter in the gas density profiles be- 
comes significantly larger in the inner regio n, which was al- 
ready noted in the previou s studies (e.g., Neuma nn & Arnaudl 
1999; Vikhlini n~etalJl999t) . 

There is also a trend for lower-temperature clusters to have 
lower gas densities and flatter profiles in the central region. 
This trend is responsible for flat gas density slopes derived for 
galaxy groups and low-mass clusters in t he previous analyses 
using the /3-model approxima t ions (e.g.. | H elsdpn <fe"Ponman 
120001 iFhioguenov et alJl200ll ISanderson et al.ll2003l) . How- 
ever, we observe that at large radii, the gas density in cold 
clusters approaches the average profile defined by T > 5 keV 
clusters. Also, the gas density profiles often steepen at r 
0.7 rsoo so that the slopes near rsoo are similar for all clusters. 
Further discussion of this issue and its impact on the hydro- 
static mass estimates is presented in AppendixlAl 

Below, we compare the concentration parameters for our 
Ptot profiles with those expected for CDM halos. Comparison 
of the gas density distributions with the results of numerical 
simulations will be presented in a future paper. 

6.1. Concentration parameters 

As discussed above, the ACDM paradigm makes a firm the- 
oretical prediction for concentrations of the dark matter halos. 
It is interesting to compare our measurements with these pre- 
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Fig. 17. — Scaled density profiles. Total density profiles are plotted within the radial range covered by the temperature profile. Gas d ensity p rofil es are extended 
to rjei (see Tablel2l. Thick yellow line shows the NFW model with C500 = 3, a typical value for CDM halos in our mass range (jj |6. II see Fig. l 1 81 . 



dictions. We define the concentration as C500 = rsoo/r s be- 
cause our mass measurements typically extend to ~ rsoo- The 
scale radius r s was determined by fitting the NFW model to 
values of the total density at six radii equally log-spaced in 
the range 0.05 r^oo - r^. The range r < 0.05 r 5 oo is excluded 
because at these radii there are separate mass components as- 
sociated with the stellar material in cD galaxies (see above); 
such components are found in a similar radial range in nu- 
merical simulat ions, which i nclude cooling and star forma- 
tion (e.g., Gnedin et al. 2004). The uncertainties for C500 are 
derived from the Monte-Carlo simulations described in §[3] 
The resulting C500 determinations are reported in Tableland 
plotted as a function of measured M500 in Figll8l 

The measurements are compared with the ex pected c(M) 
relation suggested by the cluster simulations of Dolag et al. 
(2004a, see their eq. 13 and Table 2) in the "concordance" 
ACDM cosmology, Q = 1 - £2 A = 0.3, <x 8 = 0.9 (solid 
line in Fig.1181. Note that we converted the concentra- 
tions and masses from the definitions used by Dola^gt^D 
to csn o and Mgoo used in our analysis (e.g., Hu & Kravtsov 
2003). We also show (dotted lines) the 2-cr scatter of con- 
centra tions typically found in numerical simulations, <T| n r = 
0.22 iTind l200Cjt iBullock et alJ l2Q0H IWechsler et alJ 12002 
Tasit siomi et all2004HDolag et all2004al) . 

Clearly, both the typical values and scatter of concentra- 
tions determined for our clusters are in general agreement 
with the simulation results. It can be argued that for mas- 
sive clusters, most of our measurements are slightly higher 
than the theoretical average. If this effect is real, it can be 
caused by several factors. First, our sample contains only a 



highly relaxed sub-population of galaxy clusters which are 
expecte d to sample the high tail of the concentration distri- 
bution ( Wechsl er et alJl20 02). Second, radiative cooling of 
baryons and the associated galaxy formation are expected to 
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Fig. 18. — Concentration parameters of the NFW model, C500 = '"joo/'V, as 
a function of cluster mass. Solid points with error bars are our measurements. 
Solid lin e shows the av erage concentration of CDM halos from simulations 
by Dolag et al. 1 2004a). Dotted lines show 2-<x scatter of log-normal distri- 
bution of concentrations at a fixed mass, found in simulations. 
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Fig. 19. — Total mass within rsoo (red) and r25oo (blue), as a function of X-ray spectroscopic temperature, T S n CC , and gas mass- weighted temperature, T mg . Solid 
lines show the best fit power laws with parameters listed in TableRI Dotted lines in the T spl:c plot show the best fit relations from Arnaud et al. 1 2005). Note that 
Mjoo measurements are not plotted for clusters which required large extrapolations of the temperature profiles (A262, RXJ 1 159+5531, USGC S152). 



modify the density profiles of parent CDM halos. These pro- 
cesses are expected to steepen the mass distribution in the 
inner regio ns of clusters, due t o buildup of the central clus- 
ter galaxy (Gnedinetal. 2004) and make the mass distribu- 
tion co nsiderably more spherical on scales as large as r^w 
jKazantzidis et alJl200 4 1 ). The amplitude of this effect, de- 
termined from a sample of simulated clusters described in 
iRravtsov. Naeai & Vikhhninli2005l) with two additional mas- 
sive Coma-size objects, is AC500 = 0.5 - 1 (Rravtsov et al., in 
preparation). This is similar to the possible systematic en- 
hancement of concentration we observe for M > 2 x 10 14 M Q 
clusters. 

The main conclusion from this analysis is that there is good 
overall agreement between theoretical expectations and the 
measured concentration parameter s in our cluster sample. A 
similar conclusion was reached bv lPointecouteau et alJd200"5^ 
from the XMM-Newton analysis of 10 clusters. 

7. MASS-TEMPERATURE RELATION 

A tight relation between the cluster temperature and to- 
tal mass is expected on theoretical grounds, which is in- 
deed observed at least for the hydrostatic mass est imates 
(Nevalainen et al. 2000; Horner et al. 1999; Finoguenov et al. 
120011 iXu et alJl200U iSanderson et alJl2003l) . Recent deter- 
minations of the M - T relation from Chandra an d XMM- 
New ton observations wer e presented in lAJlgn_etalJ J2001bl) 
and lArnaud et alJ d2005l) . respectively ( Al len et aT| deri ved 
masses for the critical overdensity 2500 and Arna ud et alJ also 
presented measurements at larger radii, including extrapola- 
tion to r2oo)- We can significantly improve over these previ- 
ous measurements because our temperature profiles extend to 
large radii and therefore we do not use simplifying assump- 
tions such as that T(r) is constant or polytropic. Also, we 
have high-quality X-ray surface brightness measurements for 
all clusters at radii well beyond rsoo, and do not use the com- 
mon /3-model approximation to derive the gas densities. 

We measure total masses for two often-used critical over- 
density levels, A = 2500, and A = 500. The A = 500 
level is particularly useful because it approximately delin- 
eates the inner cluster region where the bulk ICM veloci- 



ties are small and there fore the h ydrostatic mass estimates 
are meaningful (Evrard, Metzler & Navarro 1996). The cor- 
responding radi us, rsnn, is 0.5-0.67 of the virial radius, de- 
pending on Q„, (Eke et alJl996trBrvan & NormarJl998l) . The 
A = 2500 overdensity level encompasses the bright central re- 
gion where X-ray temperature profile measureme nts are fea- 
sible with Chandra even in high-redshift clusters dAllen et all 
2004). 

We do not consider the masses for lower overdensities be- 
cause this would require extrapolation far beyond the radius 
covered by the Chandra data, and because the ICM is not 
expected to be fully in hydrostatic equilibrium at large radii. 
Mass measurements are sometimes extrapolated from the in- 
ner region to, e.g., r2oo assuming the NFW model for the 
matter density profile. Such extrapolations are highly model- 
dependent and lead to underestimated measurement uncer- 
tainties for M2oo- It is more appropriate to scale the theoreti- 
cal models to rsoo, where direct measurements are now avail- 
able. Scaling of mass of CD M halos to a ny overdensity level 
is straightforward (Hu & Kravtsov 2003). 

The measured values of M2500 an d ^500 are reported in Ta- 
ble[3] We do not compute M500 for three clusters for which the 
temperature profile measurements have to be extrapolated too 
far (A262, RXJ 1 159+5531, and USGC S152). For A2390, 
M2500 should be treated with caution because ^500 is near the 
boundary of the central non-hydrostatic region (see §@J. 

In Fig^] we plot the measured masses as a function of 
the X-ray spectroscopic and gas mass-weighted temperatures. 
The self-similar evolution of the normalization of the M - T 
relation is expected to follow 

M A /T 3/2 cc E(zT\ E(z) = H(z)/H , (11) 

(e.g. JBrvan & Normanll9 98). Even though the effect of evo- 
lution is small within our redshift interval, we applied the cor- 
responding corrections by multiplying the measured masses 
by E(z) computed for our preferred cosmology, thus "adjust- 
ing" all measurements to z — 0. 
We fit the observed M - T relation with a power law, 

M = M 5 (775 keV) a . (12) 
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TABLE 4 

Power Law Fit to Mass-Temperature Relation 













T 

1 mg 




Overdensity 


m 5 , r'M Q 


rs, ft -1 Mpc 


a 


M 5 , h- l M e 


rs, h Mpc 


a 


A = 500 

A = 2500 


(2.93 + 0.16) XlO 14 
(1.28 + 0.05) x 10 14 


0.796 ± 0.015 
0.354 ± 0.005 


1.61 ±0.11 
1.64 ±0.06 


(3.41 ±0.18) x 10 14 
(1.48 ± 0.07) x 1() 14 


0.837 ±0.015 
0.371 ±0.006 


1.51 ±0.10 
1.58 ±0.07 



Note. — We use a power law fit of the f orm h (z) M = Ms X (775 keV) a , where M = M500 or M2500 and temperatures are eithe r X-ray 
spectral or gas mass-weighted averages (see § 13.51 . Scaling of the corresponding overdensity radii with temperature is of the form 1131 . 



The relation is normalized to T — 5 keV because this is ap- 
proximately the median temperature for our sample and there- 
fore the estimates for M5 and a should be nearly independent. 
The fit is performed usin g the bisector modification of the 
lAkritas & Bershadvl i 19961 and references therein) linear re- 
gression algorithm that allows for intrinsic scatter and nonuni- 
form measurement errors in both variables. The uncertain- 
ties w ere evaluated by bootstrap resampling (e.g., Press et aTl 
1992), while simultaneously adding random measurement er- 
rors to M and T. 

The best fit slopes and normalizations for different variants 
of the M - T relations are reported in Tabled We find flatter 
slopes (a - 1.5 - 1.6) than many previous studies (typically, 
ot = 1.7 — 1.8 when low-temperature clusters were included). 
Some of the difference can be traced to slightly different def- 
initions of the cluster temperature in individual studies and 
different procedures for scaling the measurements to z = 0. 
However, the major effects are accurate measurements of the 
temperature gradient at large radii and correct modeling of the 
steepening in the gas density profiles at large radii. The com- 
parison with other works discussed in detail in AppendixlAl 

The M - T relation implies the following scaling of the 
overdensity radii with temperature 



r A hE(z) = r 5 (775keV) 



a/3 



(13) 



with coefficients rg provided in Tabled 

As seen in Fig^| the scatter of the individual M and T 
measurements around the best fit power law approximations 
is very small. Note that in our case, unlike many previous 
studies, this is not just a trivial consequence of the approxi- 
mate similarity of the cluster X-ray surface brightness profiles 
(N eumann & Arnaud 1999) because we do not use overly- 
constrained models for T(r). To characterize the scatter, we 
compute the rms deviations in mass and subtract the expected 
contribution from the measurement errors, 



(6M/M) = 



1 



N-2 



Z 



(Mi- M 5 (Ti/5) a r - AM, 



(14) 



where AM, are measurement errors and SM/M is the esti- 
mated scatter in the relation (the l/(N - 2) factor accounts 
for 2 degrees of freedom in the power law fit). We find that 
the observed scatter is consistent with zero, and the 90% up- 
per limits for intrinsic scatter are 6M/M ^ 0. 15 for both M500 
and M2500, and for both definitions of the mean temperature, 
r spec and T mg . Note the scatter is likely to be larger for the 
whole cluster population including non-relaxed objects. 

Our normalizations of the M - T relation are higher than 
most of the previous X-ray determinations based on ASCA 
and ROSAT analyses (see Appendix lAl for detailed discus- 
sion). There is, however, a very goo d agreement with the 
XMM-Newton measurements by Arnaud et al. (2005, their re- 
sults are shown by dotted lines in Fig|19ll. although this ap- 



pears to be a result of a chance cancellati on of systematic dif- 
ferences between our analyses (Appendix lA.2> . 

Our normalization of the M - T relation is also in good 
agreement with those deri ved from r ecent high-resolution nu- 
merical simulations (e.g., Bor gani et alJl2004h that attempt to 
model non-gravitational processes in the ICM (radiative cool- 
ing, star formation, feedback from SN). Detailed comparison 
of our measurements with the results of numerical simulations 
will be presented elsewhere. 

The primary goal of observationally calibrating the M - T 
relation is for use in fitting cosmological models to the clus- 
ter temperature function. Given the good agre ement of our 
result s with some other recent measurements jArnaud et alJ 
120051) and results of realistic cluster numerical simulations, 
it is tempting to conclude that observational determinations 
of the M - T normalization have finally converged to the 
true value. We, however, caution against direct application of 
our normalizations to fitting the published cluster temperature 
functions f M arkevitchll 19981 lHenrvll200fl: llkebe et al]l2002l) 
for several reasons. First, the definitions of the mean cluster 
temperature in these papers is not identical to ours. Second, 
our determination of the M — T relation uses only the most re- 
laxed clusters at the p resent epoch and it can be significantly 
different for mergers (Randall et al. 2002). Even for relaxed 
clusters in the present sample, our analysis neglects poten- 
tial deviations from hydrostatic equilibrium (e.g., caused by 
ICM turbulence) or deviations from spherical symmetry. A 
detailed study of these effects will be presented in Nagai et 
al. (in preparation); they can lead to ~ 10% underestimation 
of the total masses in our analysis ("e.g.. [ Ev rard et a lJl!996t 
iRasia et alJ2004tlKav et alJ2004l: lFaltenbac her et all2005l) . 

8. GAS FRACTIONS 

Finally, we present determinations of the ICM mass frac- 
tions in our clusters. Such measurements, especially within a 
large fraction of the cluster virial radius, are of great cosmo- 
logical significance because they provide an independent test 
for Q.m assuming tha t the baryon fractio n in clusters is close 
to the cosmic mean (White et al. 1993b). 

The derived ICM mass fractions, / gas , as a function of radius 
are shown for individual clusters in Fig.l3*t-[l5b. We present 
results for both enclosed gas fractions (continuous lines) and 
local fractions within several spherical shells (crosses). In al- 
most all clusters there is a significant increase of / gas with 
radius, in qualitative agreement with a number of previ- 
ous measu rements that used spatial l y-resolved temperatures 
(e.g.. IMarkevitch & yikhlininlll997l: IMarkevitch etaflll999l: 
iPratt & Arnaudl2002UAllen et all2004l) . 

Comparison of observed gas fractions for individual clus- 
ters as a function of overdensity (Figl20> gives further in- 
sights. At A < 10 4 , the gas fraction increases with radius as a 
power law of overdensity. Except for maybe in three clusters, 
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no flattening of / gas (r) is observed, at least within r5oo- The 
average ratio of / gas within r 2 500 and r 50 o is /2500//500 = 0.84 
in massive clusters. 

There is also a significant trend of increasing / gas with clus- 
ter mass. The effect is the strongest at small radii, but even 
within 7-500, / gas (M) is not constant. In Fig|^ we plot the 
enclosed gas fractions within r25oo and rsoo as a function of 
cluster temperature. The gas fraction at r25oo increases ap- 
proximately linearly with temperature from /2500 — 0.04 for 
three T < 2 keV clusters to /2500 - 0.10 - 0.1 1 for the most 
massive clusters. Another possibility is that there is a flatten- 
ing at /2500 - 0.09 at T > 5 keV and that the measurement 
for the highest-temperature cluster, A2390, is strongly biased 
(see §0}. We cannot distinguish these possibilities because of 
significant object-to-object scatter and the small sample size. 

The global baryon fraction in the Universe is constrained 
by CMB observations to be Q/,/Q m = 0.175 + 0.023 
jReadhead et alJ 120041 ISpergel et alJEOOl . Therefore, the 
observed gas fraction within f25oo> even m the most massive 
of our clusters, is significantly lower, by a factor of ~ 0.6, 
than the cosmic mea n. A similar le vel of baryonic deficit was 
previously noted by Ettori (2003) and its syste matic varia- 
tion w ith mass was obtained i ndirectly by lArnaud &"E vrard 
(1999) and iMohr etafl (H999). This deficit, and cluster-to- 
cluster variations of /2500 can, at least in part, be explained by 
conversion of ICM into stars. Contribution of the stars to the 
total baryon budget should be most important in the cluster 
centers because of the presence of cD galaxies. Indeed, the 
largest cD galaxies have K-band luminosities (1 -2)xlO 12 L 
dLin & Mohrl2004l). Assuming a st ellar mass-to-light ratio in 
the K-band of ~ 1 dBell et alJ20 03). we can estimate that just 
the cD, not counting other galaxies, can contribute from 0.07 
to 0.01 to the baryon fraction within r25oo, depending on the 
cluster mass. Therefore, the stellar mass in the cluster cen- 
ters is significant and it should be determined individually in 
each cluster. Another important process that can affect / gas in 
the cluster centers — and otherwise break the self-similarity 
of the gas propertie s — is energy output from central AGNs 
JNulsen et all2005l) . 

The stellar contribution to the total baryon mass and rela- 
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Fig. 20. — Enclosed gas fraction as a function of overdensity, denned with 
respect to the critical density at the cluster redshift. Dashed lines show the 
range Q.t / fi„, = 0. 1 75 ± 0.023 constrained by CMB observations. Points with 
error bars show typical measurement uncertainties at two radii. 



I I I I I 1 I 


CMB ] 






: I 


} + 4 + + + ; 


: f * 


• 


♦+ 







2 4 6 8 10 

■^spec > keV 



Fig. 21. — Observed gas fractions within r25oo and rsoo (solid and open 
circles, respectively). We do not compute /500 f° r A262, RXJ 1159+5531, 
and USGC S 152 because this requires large extrapolations of the mass model. 
The / 250 o value for A2390 (T = 8.9 keV), can be biased (§@. 

tive energetics of non-gravitational processes should be less 
important within larger radii. Indeed, we observe that / gas in- 
creases between r25oo and rsoo, by factors of 1.2-1.4. Also, the 
trend of /500 with temperature is weaker than that for /isoo- It 
is consistent with both a linear increase of / gas with T and 
flattening at T > 5 keV. 

Note, however, that typically, 40-50% of the total mass 
within rsoo is contained within the region r < r 2 5oo affected by 
stellar contribution and non-gravitational heating by AGNs. 
One can hope that a greater uniformity is observed if this cen- 
tral region is removed entirely from the calculation of / gas . 
In Fig|22j we plot gas fractions measured in the r25oo-'"500 
shell. The gas fractions in the shell increase still further rel- 
ative to /5oo- In fact, they become consistent (albeit within 
larger uncertainties) with an average of / gas =s 0.13 for all 
clusters, except for the lowest-temperature object (MKW 4). 
A reasonable correction for the stellar co ntributio n, 10- 
15% of the ICM mass (e g iLin. Mohr & Stanford l2003t 
ISanderson & Ponmanl 12001 IVoevodkin & Vikhlininl 120041) 
should then bring the gas fractions within the shell r 2 soo- 
r 5 oo into agreement with the Universal value determined from 
CMB studies. 

To summarize, the observed gas fraction shows smaller 
variations between individual objects, and is closer to the Uni- 
versal value, when determined at larger radii. This is in line 
with the general tendency for our clusters to become more 
self-similar at large radii, as manifested by the 3-dimensional 
temperature and density profiles, and the M - T relation. The 
regularity of / gas determined in the shell ^500-7500 gives hope 
that these measurements can be used to determine Q m using 
the classical baryon fraction test ( White et alJl993bl) . We de- 
fer application of this test to a future work because this in- 
volves small but important corrections, such as stellar contri- 
bution, baryon depletion, accuracy of hydrostatic assumption 
etc., which are beyond the scope of this paper. 

We finally mention that gas fractions within r 2 soo derived 
from Chandra observations of relaxed clusters were recently 
used for cosmo logical constraints in a series of papers by 
lAllen et alJ ( 12004 and refer ences therein). There are signifi- 
cant differences in the /2500 derived in our analysis and those 



CLUSTER MASS PROFILES 



15 



-i — i — i — i — | — i — i — ' i 



0.20 



0.15 



0.10 



CMB 



0.05 - 



0.00 



10 



Fig. 22. — Observed gas fractions within the shell rqsgo— rsog. We do not 
compute / ga8 for A262, RXJ 1159+5531, and USGC S152 because this re- 
quires large extrapolations of the mass model. 

reported bv lAllen et all as outlined in AnnendixlBl 

9. SUMMARY AND FUTURE WORK 

We presented total mass and gas profiles of low-redshift, 
relaxed clusters using the best available observations with 
Chandra. The cluster sample (13 objects) spans a range 
of temperatures of 0.7-9 keV and masses M500 = (0.5 - 
10) x 1O 14 M . The total masses are derived assuming hy- 
drostatic equilibrium of the ICM in the cluster potential. Our 
modeling method makes few additional assumptions and al- 
lows us to reconstruct the cluster properties relatively model- 
independently, thus avoiding biases and producing realistic 
uncertainties on the derived quantities. The main results are 
summarized below. 

1. The shape of the total density profiles in our clusters 
and their scaling with mass are in good agreement with pre- 
dictions of ACDM model. They approximately follow the 
universal density profile with the concentration expected for 
CDM-dominated halos in this cosmology. The gas density 
and temperature profiles at large radii are also nearly self- 
similar, although in the inner region (within « 0.1 r 5 oo) there 
is significant scatter and systematic trend with cluster mass. 



2. Correspondingly, we find that the slope of the mass- 
temperature relation is in good agreement with the simple 
self-similar behavior, M500 K T a where <x= (1.5— 1.6) + 0.1, 
if the average temperature is measured at radii not affected by 
the central cool core. We derive an accurate normalization of 
the M - T relation for relaxed clusters. Our normalization is 
~ 30% higher than most previous X-ray determinations. 

3. The gas density profiles are generally not described by 
a yS-model or its common modifications, even at large radii. 
The profiles steepen continuously and this behavior can be ap- 
proximated as a smooth break in the power law index around 
(0.5-0.7) rsoo. Near rsoo, the effective slope of the gas density 
profile, — I d\ogp g /d\ogr 0.78, with no detectable depen- 
dence on cluster mass. 

This behavior of the gas density at large radii is missed by 
/3-model fits because they are primarily sensitive to the data in 
the bright inner region. Insufficiently accurate modeling of the 
gas distribution at large radii, in addition to using polytropic 
approximations to T(f), explains the lower normalizations of 
the M - T relation derived in previous studies. 

4. We present accurate measurements of the gas mass frac- 
tion as a function of radius. We observe strong systematic 
variations of f gas both with radius and with cluster mass. The 
gas fractions within r25oo are significantly lower than the Uni- 
versal baryon fraction suggested by the CMB observations. 
However, the trends become weaker and the absolute values 
°f /gas are closer to the Universal value at r > r^soo- 

In future work, we will use these accurate measurements 
of the gas fractions to constrain Q m . This requires the in- 
clusion of small but important effects, such as stellar mass, 
baryon depletion, and correction for biases in the mass mea- 
surements, which are beyond of the scope of this paper. We 
will also present a detailed comparison of the Chandra results 
on the cluster mass, temperature, and / gas profiles with high- 
resolution cosmological simulations. 
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APPENDIX 

A. COMPARISON WITH PREVIOUS DETERMINATIONS OF THE M-T RELATION 

The main differences of our analysis with the previous work on cluster masses is in the more accurate measurements of the 
gas temperature and density gradients at large radii. Comparison between different studies is facilitated by formulation of the the 
hydrostatic equlibrium equation (Q) using effective slopes of the density and temperature profiles. 

Let us define the effective bas density slope as /3 e ff = —h dlogp/dlog r; for the y8-model, y6 e ff = /3 (r/r c ) 2 /(l + (r/r c ) 2 ) — > /? at 
large radii. The equivalent quantity for the temperature profile is ft, = -1/3 c/log 77 of log r. For polytropic parameterization of 
the temperature profile, y8 e ff and f3 t are related via 

-77-^— = (r - i)~77~^ or A = (y-DAff (Al) 

a log r a log r 

The hydrostatic equilibrium equation (7) can now be rewritten as M(r) oc r x T(r) x (J3 e ff + /3,). To compute the overdensity 
mass, we solve equation of the type M(r&)/r 3 A = C, or r A 2 T(r) (j8 e ff + f3 t ) - C. Therefore, the overdensity mass estimate scales as 

M A cc r 3/ 2 (T(r)/T ? /2 (J3 eS + Pd V \ (A2) 
where To is an average temperature, and normalization of the M-T relation scales as 

A = M A /T 3/2 cc (T(r)/T ) V2 (J3 e « + #) 3/2 (A3) 
The quantities /3 e ff and f3, derived for our clusters at rsoo are shown in Fig lAll 

A. 1 . ASCA Measurements by Nevalainen et al. and Finoguenov et al. 

iNevalainen et alJ i2000l) andlFinoguen ov et alJ (|2001 ) use ASCA temperature profiles which agree with our Chandra measure- 
ments. However, their normalizations for the Msm — T re lation for T — 5 keV clusters are lower — 2.2 x lO 14 /^ 1 M and 
2.3 x 10 14 h^ 1 M Q in Nevalai nen et alJ and Finogueno v et all respectively, compared with our value of 2.93 x 10 14 h~ l M G (Ta- 
ble |4] for r spec ). There are subtle differences between these studies in how the average temperatures are defined and how the 
measurements are scaled to z = 0. The most significant effect, however, is steepening of the gas density profile at large radii 
that is clearly present in high -quality C ha ndra data and was hard to detect with earlier X-ray telescopes. This steepening was 
effectively missed by Nevalai nen et alJ and Finogueno v et alJ because they used pur e jS-model fits determined mainly by the data 
in the inner region s. Indeed, the average /3 e ff at rsoo is 0.78 for our clusters (Fig. lAH . while i t is » 0.65 f or T > 3 keV clusters in 
Finog uenov et alJ (see their Fig. 5). The average /3, for our clu sters is 0.17, while the Finog uenov et alJ value is 0.1 1 from their 
average polytropic index of y » 1.17 for hot clusters (see eq. IA1I ). These differences in the slopes should lead (eq. IA2I 1 to a 
mismatch in mass estimates by a factor of (0.78 + 0.15) 3 ^ 2 /(0.65 + O.ll) 3 ^ 2 1.35, explaining the offset betwe en our M-T 
normalizations. We note here that undetected steepenin g of the gas density profiles at large radii was suggested by Borgani et al. 
(2004) as a possible reason for low normalization of the Finogue nov et alJ M - T relation. 
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We find no detectable trends in either y6 e fr or /3, with cluster temperature (Fig. lAH . The gas density profiles of cool clusters are 
indeed flatter in the inner region but they steepen significantly at r ~ 0.5 r^oo (Fig.ll7>. T his a bsence of trends in y8 e ff and /3, is in 
fact the main reason why our slopes of mass-temperature relation are close to 1.5 (cf. eq. IIA2I ). In contrast, earlier studies based 
on y6-model fits to the X-ray brightness p rofiles consistently found very flat density slopes for 1-2 keV clusters. For example, the 
average /3 e g for such clusters is ^ 0.5 in Finogue nov et all while for T ~ 10 keV clusters, they find /3 e g w 0.7. It follows from 
eq. IA2I that such a trend in /3 e g should steepen the slope of the M - T relation from a = 1 .5 to a ~ 1 .72. 

A.2. XMM-Newton Measurements of Arnaud et al. 

Our normalization of the M - T relation is very close to that derived bv lArnaud et alJ J2005I) using XMM-Newton data. Given 
the systematic difference in the temperature profiles at large radii between our two studies (see the discussion in §[5}, this is not 
expected and merits a clari fication. In the Arnaud et al. works, temperature profiles are generally interpreted to be consistent 
with constant at large radii ( Pointecouteau et al. 2005). Therefore, the quantities to use in eqs. (A2-A3), /3, - and T(r)/To 1 
near r^ for the XMM-Newton masses. E. Pointecouteau and M. Arnaud kindly provided the average gas density slope for their 
sample, (y6 e fr) = 0.71, only slightly below our value. For our sample, the relevant qua ntities at r^p p are (J3 e s) = 0.78, (J3 t ) = 0.17, 
and (T(r)/To) = 0.67, where Tq is the spectroscopic average temperature (see § I3.5> . From eq. IA3I we would then expect the 
XMM-Newton M500 - T normalization to be a factor of 1.18 higher than our value, while in fact it is slightly lower (Figll9t. 

The reason is that for the M -T relation, Arnaud et al. used cluster masses derived by fitting an NFW model to the data within 
~ '"looo. their maximum radius of observation. This gives systematically lower masses at r > rsoo than the values obtained by 
direct hydrostatic derivation using extrapolation of the isothermal temperature profiles (E. Pointecouteau & M. Arnaud, private 
communication). In effect, the masses from the NFW fit imply declining temperature profiles at large radii, such as those observed 
by Chandra. To summarize, the agreement between our M — T normalizations is somewhat a coincidence. 

Note also that th e XMM-Newton cluster temperatures are systematically lower than those from Chandra (see e.g., 
Vikhl irun et al.1120051) '). but this does not change the normalization of the M - T relation as it moves the clusters along the 
relation. 

We also used the results from the iPointecouteau et alJ (120 05) sample to check our mass derivation algorithm. The published 
XMM-Newton temper ature and density pro files were used as input to our procedure. The obtained mass profiles were nearly 
identical to those from Pointecouteau et al. 



B. COMPARISON WITH / GAS MEASUREMENTS OF ALLEN ET AL. 

Gas fractions within ^2500 d erived from Chand ra observations of relaxed clusters were recently used for cosmological con- 
straints in a series of pap ers bv | Allen et alJ (12004. and references therein). There are significant differences in the / gas ,2500 derived 
in our analysis and in Allen et al., as outlined below. 

For all our T > 5 keV clusters except A2390, we derive l ower values of / gas . Our average value for these clusters is / gas = 
0.091 + 0.002, compared to an average of 0.117 + 0.002 in Allen et al. (same temperature range and same cosmology). The 
same ~ 25% di fference holds for the four clusters common in both samples, A2029, A478, A1413, and A383 (see Table 2 in 
lAllen etall2004l and our Table QJ. 

We also do not observe flattening of f ^jr) at radii withi n 7-2500 as reported, e.g., in lAllen. Schmidt & Fabiar] 12002). Note, 
however, that a larger sample presented in Allen et al. ( 2004) shows more variety in the behavior of f gas (r) profiles. 

Profiles for individual clusters in the I Allen et alJ sample were either not published or published prior to significant Chandra 
calibration updates. Therefore, we are unable to perform a more detailed object-to-object comparison. 



